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ABSTRACT 

Constrained effective potentials in hot gauge theory give the probability that a con- 
figuration p of the order parameter (Polyakov loop) occurs. They are important in the 
analysis of surface effects and bubble formation in the plasma. The vector potential ap- 
pears non-linearly in the loop; in weak coupling the linear term gives rise to the traditional 
free energy graphs. But the non-linear terms generate insertions of the constrained modes 
into the free energy graphs, through renormalisations of the Polyakov loop. These inser- 
tions are gauge dependent and are necessary to cancel the gauge dependence of the free 
energy graphs. The latter is shown, through the BRST identities, to have again the form 
of constrained mode insertions. It also follows, that absolute minima of the potential are 
at the centergroup values of the loop. We evaluate the two- loop contributions for SU(N) 
gauge theories, with and without quarks, for the full domain of the N-l variables. 
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1. Introduction 



Study of high-temperature QCD 1 ) has become urgent because of the building of the 
relativistic ion collider (RHIC). Apart from this, there is the vast area of early cosmology, 
where hot QCD has come into play. At these temperatures the Z(N) symmetry of quarkless 
QCD is broken spontaneously leaving us with N ordered phases. Adding quarks breaks 
this symmetry explicitely, but the deconfmement transition is still there. 

Recently there has been interest 2 ) 3 ) in the surface tension that occurs in QCD with 
quarks only as an order-disorder interface, and, subsequently, in the surface tension be- 
tween the different ordered phases in quarkless QCD. These ordered phases are charac- 
terized by the order parameter (the Polyakov loop) having a Z(N) center group value. 
Numerical work on this preceded our analytic approach 4 ) 5 ), and both are in reasonable 
agreement. 

One key ingredient in these calculations is the profile p that the Polyakov loop p 
develops in between two Z(N) phases. The other key ingredient is the probability that a 
given profile occurs. This quantity is called a "constrained" effective potential. 

In this paper we will calculate the constrained effective potential U, that is associated 
with the probability that a given constant mode p = y f dx p(x) appears : 



exp[-V0U(t)] 



fDA6(t-p)ex.p-j*S(A) 
fDAexp-^S(A) 



U(t) as a function of t is computable in perturbation theory. It is accessible to Montecarlo 
simulations and is gauge independent. We will give the perturbative rules for computing 
it ; these rules include the ad-hoc prescription of ref 12. 

In a previous paper 5 ) we did compute the constrained effective action for a special one 
parameter set of the Polyakov loop, in the absence of quarks. This special set is called the 
" q valley" . Its physical significance is due to the fact that it vehicles the tunneling effects 
giving rise to the surface tension between two ordered vacua 5 - 1 . 

In this paper we will give the full result, including all values of the Polyakov loop 
outside the q valley, and including quarks. We show (in section 3) that the ad hoc insertion 
method in 12 ) and 5 ) follows from starting from the constrained path integral. Our results 
are useful for calculations of transition probabilities between various stable and unstable 
minima that appear in the quark case, e.g. for the case of three colours and six quarks 15 ). 

Apart from this we have been motivated by two reasons 
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i) To lift some of the uncertainties in the literature on the pure gluon results. Those 
readers that are only interested in the end results can consult eqns 4.2, 4.3 (pure glue) 
5.2 and 5.4 (fermions). Our pure glue results are in agreement with ref 5 eqn 5.1, 5.13 
and 5.14. 

ii) To gain some insight into the gauge-cancellations using BRST identities generalized 
to a non-trivial Polyakov loop background, and the accompanying stabilisation of the 
Z(N) minima due to extra vertices induced by the background. There is some overlap 
with ref 13, but our methods are complementary. The gauge artifact cancellations are 
shown in 2.35 and 3.18. 

As mentioned, there has been related work by various authors, with different motivations 6 ) 7 ),! 
and different outcomes. With our method we compute in various gauge fixing choices and 
find the same result. This fortifies our confidence in our results. Besides, we find that the 
Z(N) vacua stay stable to two loop order ; if quarks are present the stable vacuum stays 
at where the Polyakov loop takes the value one, where as metastable minima remain at 
roughly the location of the non-trivial center group values of the loop (see section 5 and 
6). 

In section 2 we explain in some detail the method of calculation with as main result 
eq. 2.35. The expectation value of the Polyakov loop to one loop order is computed in 
section 3, resulting in eq. 3.17 and 3.18. In section 4 the gluonic 2-loop results are given in 
4.2 and 4.3.. In section 5 the fermionic results are analysed. Section 6 contains conclusions 
and outlook. A few appendices elaborate on the results. 
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2. The method to compute the effective action U. 



This section is rather long so it is divided into six subsections. Section 2. a defines the 
effective action, section 2.b gives the saddle point method for the perturbative evaluation 
and section 2.c deals with the one loop results. In section 2.d the two loop results are 
explained; it is here that the intricacies of the non linearities of the Polyakov loop enter. 
How the final answer is gauge choice independent is detailed in sections 2.e and 2.f. 

2. a. Effective potential U. 

We will work in a finite box V = L s in a heat both of temperature T = ; in the box 
we have an SU(N) Yang Mills field, and we work on the usual Euclidean formulation with 
periodic b.c.'s in the time direction. This means that the unitary NxN matrix P(Aq(x)) 
defined by 

rP 

P(A (x)) = V(expi / cItA (xt)) = lim j [ expiATkA (x, r^) 2.1 

J0 k=l 

- V being the path ordering, rt+i = tu + At, At = (3/n - is transformed by a periodic 
gauge transformation O like QP(A )Q~ 1 , and has therefore gauge invariant eigenvalues. 
Another way of saying this that : 

hiAoix)) ^^Tr P(AoW),-~,tN-i(M*)) ^ ^ Tr P N ~\MS)) 2.1(a) 

are gauge invariant. For Tr P(Aq(x)), the "Polyakov loop" we have a straightforward 
interpretation once it is averaged in the Euclidean path integral formalism : 

^(M*))) = JZMexp _j^ (j4) 2-2 

This average is the exponential of the free energy of an infinitely heavy quark (i.e. a 
fermion in the fundamental representation of SU(N)) at the point x. {ti(x)) depends on 
the spatial boundary conditions. This simple relation between the average and the free 
energy F q of a single quark reads: 

(t 1 (A (x))) = exp-PF q 2.3 

From this relation ^ follows the use of < ti(x) > as an order parameter. 
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Averages of traces of higher powers of the loop 2.1 t k (x) can be gotten by suitable 
linear combinations of free energies of "quarks" in higher representations of SU(N). Taking 
the SU(3) case as an example, one can relate 

< TrP(A ) 2 >=< TrP 6 (A ) >-< TrP(A ) > 

with P 6 the Polyakov loop in the sextet representation of SU(3). From now on we will be 
interested in spatial averages tk = y f dx tk(A (x)), and we will define, following Fukuda 
et al 8 ) 9 ), the constrained effective potential U by: 

r N_1 1 

e W -/3V(U(t 1 ...t N . 1 )+F) = / DA ][ S(t k - t k ) exp - — S(A) 2.4 

J k=i g 

The left hand side is up to normalisation the probability that a given set of fixed numerical 
values t\ . . . tjv-i f° r the corresponding spatial averages appears in our system 20 . The 
normalisation is in terms of the free energy F, which one gets by doing the path integral 
2.4 without the constraints. In ref 9 it was argued that this effective potential is numerically 
a very useful quantity, and easier accessible than say the usual free energy in the presence 
of sources J\ . . . Jjv-i for the tk '■ 

1 N ~ x 
exp -f3VF(J t . . . t/jY— i ) = / DA exp --S(A) + V ^ J k ■ t k 

J g k=i 

This free energy is obviously related to 2.4 through the Laplace transform : 

exp — 13V F( Ji . . . Jjv-i) = J dt\ . . . dtN-i exp V ' J\t\ + . . . + V J^-itM-i 

exp-pVUfa.-.tN-!) 2.5 

U and F are clearly gauge invariant. 

Consider the infinite volume limit in eqn 2.5. The ensuing saddle point equations give 
U the same dependence on the t k as the effective potential G = F — §j.J has on the 
expectation values of the tk- In fact one has 

G = U + F(0) 2.6 

in the infinite volume limit. The reader may ask at this point: why not discuss the effective 
potential instead of the constrained effective potential? A first reason for preferring the 
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latter is that it is given directly in terms of a path integral, eqn 2.4. Second, in a finite 
volume there are differences, which may become significant in perturbation theory, when 
there are infrared divergencies. Third, the study of surface effects necessitates the use of 
the constrained effective potential. 

Some authors 13 ) prefer to discuss the eigenvalues of the loop, rather than the traces. 
The Jacobian between the two is a van der Monde determinant, as explained in appendix 
D. We prefer traces, because they are easier to handle when taking higher order effects 
into account (see section 3). 

The next subsection deals with the perturbative evaluation of the effective potential. 

2.b. Perturbative evaluation. 

As in ref 5) we will Fourier transform the ^-constraints. This will give us N — 1 
variables A&, k = 1 . . . N — 1, with : 



JV-l „AT-1 



n / n^ e ^ (tfe_ffe)) 2 - 7 



k=l " fe=l 



In the following we will drop the index k for notational simplicity, but indicate the depen- 
dence of t on the potential Aq. Substituting 2.7 into 2.4 leaves us with a path integral in 
terms at the gauge invariant variables A and the gauge potentials A : 

exp-(3VU(t) = J dxj DAexpi\(t -i(A )) - ^S(A) 2.8 

To keep the bookkeeping of our degrees of freedom precise, we have to introduce 
boundary conditions, e.g. periodic boundary conditions. 

We will search for a saddlepoint A = 6, A = B and set therefore : 

A fl = B fl + gQ^ 
\ = b + gq 2.9 

where Q and q are quantumvariables, expand the exponent in 2.8 and require that linear 
terms cancel. The integration over the quantum variables in 2.9 will leave us with a 
functional Z(b, B) which stays invariant under a gauge transformation 

A a = 041T 1 + fidfi- 1 
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The latter can be rewritten in obvious way as 



B n = OBfi" 1 + fidfT 1 2.10 

together with 

Q n = QQQ- 1 

It is then clear that indeed 

Z(b,B) = Z(b,B Q ) 2.11 

The path integral 2.8 still needs gauge fixing ; in order to retain property 2.11, we will 
take the familiar background field gauge fixing : 

S 9 .f. = ^Tr(D,(B)Q,) 2 2.12 

with the Faddeev-Popov term : 

r]D^B)D^B + Q)iv 2.13 

Using 2.1 and 2.9 we find for the linear term in an obvious short hand notation : 

i(t - i(B))q + ibt^(B)Q - SW(B) • Q = 2.14 

This gives three saddle point equations, one for the q's (2.15(a)), one for those Q's, that 
couple to t^(B) and S^(B) (2.15(b)) and finally those Q's, that couple only to S^(B) : 

t-t(B) = 2.15(a) 
ibi^(B)- S^(B) =0 2.15(6) 
S^(B) = 2.15(c) 

Since 2.11 determines B up to a gauge transformation we choose a solution B that is 
computationally most convenient : 

S M = Co M0 2.16 

C is a space time independent diagonal matrix in the Lie algebra of SU(N) ; this satisfies 
2.15 with 6 = 0. We emphasize this is just a choice, and will come back to it in section 
2.e, where the gauge variance of our results is discussed. 
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So the eigenvalues Cj(z = 1 . . . N) of C obey a constraint 



N 



^C, = 2.17 



i=i 



and are related to the N — 1 (k = 1, . . . , N — 1) through = ^Tr expikC '(3. 

The advantage of this choice for the saddle point lies in the easy diagonalisation of 
the quadratic part of the action : 

5(2) =iq^ J dxdTTri {1 \C)Q {xT) 

+ J dxdr TrQ^-D\C) + (1 - £)D, i (C)D v (C))Q v + 

2 J dxdrTrf](-D 2 (C))u 2.18 

Dn(C) is the covariant derivative : D^{C) = + 5^q[C,. 

Let us diagonalize the last two terms in 2.18, in terms of the colour basis (A lJ ) n>m = 
■^j^in, 8jm and X d = -^diag(l, . . . , 1, 1 — d, . . . , 0), Tr = 0. Details on this Cartan 
basis are discussed in appendix D. 

The Fourier transform of the Q's is defined as : 

QM = LJd2j* dre-^-^Q^x, r) 2.19 

with 

P l 1 2 fs°}w = 0,1,2,... 2.19(a) 
We write Q^(p) in terms of colour components as : 

JV 

QM = Y.Q» xij + J E®» xd 2 - 20 

and find for in its diagonalized form, restoring the colour indices: 
N-i v 

s {2) = * E «^i(c)QS(o) + -t^ gj?(p)((p y ) v - (i - dp'frftQi'i-p) 

v N 

d=2 n 
TV 

+ vrj2J2 f]d (p)p 2ujd (-p) 2 - 21 

d=2 n„ 



In the first term the matrix tj j = Tr exp(ilC)Xd results from expanding the 1-th power 
of the Polyakov loop to first order in the quantum field Qq(0). The symbol p^ stands for 
Po = Po + Ci — Cj, and p l £ = pk- This shift in the momentum is due to the covariant 
derivatives in 2.18. 



is due to our periodic boundary conditions. We can use e.g. twisted b.c!s and overcome 
this problem to get the same thermodynamic limit as we will get here by ignoring the 
problem. 

The N — 1 constraint variables q couple only to the N — 1 zero-momentum diagonal 
Qo variables, Qq(0), in the first term of 2.21. 

2.c. One loop result 

The one loop result is obtained from the expression for in 2.21, by substituting 
it in the exponent of 2.8. 

We have first to integrate out the q variables. This will give us delta function con- 
straints on the N-l Qq(0) variables, of the type 5(Qq(0)), and the one loop result is: 



The prime on the trace means we left out the eigenvalues corresponding to these N — l 
modes. 

We are interested in the thermodynamic limit V — > oo. Then these N — 1 constraints, 
as argued in ref 5), are not important. 



Let us note that 



r) k (p = 0),uj k (p = 0) and Q^ip = 0) are non-Gaussian variables. This 



VpTjV (t) = \tt' log(4^ - (1 - Opfp)!) ~ Tr' \ogp% 



2.22 



For completeness we mention the result ^ : 




2.23 



where Cij = 



and have to be taken mod 1; the t are given in terms of the diagonal 



matrix C 2.16 : 




Note that is independent of the gauge choice, as it should. So far, all our labours did 
not lead to anything new. This changes in the next subsection. 



2.d. Two loop results 
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For the two loop results we need to expand the action and the constraint to order g 2 . 
To avoid clutter in the formulas we will stop indicating the colour degrees of freedom on 
q, t (2) (C) etc. . . . This leads to : 

s {3) ^(|t (2) (^)-Q!+^ (3) (^)-Q!) 

+ f S£l • Q 3 + f^nl • Q 4 + 2(D,(C)f))[QM 2.25 

The bars mean the average over the space- volume V. 

is standard, except for the first term linear in q. Had the constraint been linear 
in the potential, the second and third derivative of i(C) would have been absent. Then 
we would have had - as in the one-loop case - as only effect from the g-integration the 
5-function constraint on the N — 1 Qq(p = 0) variables ; hence, in the thermodynamic 
limit, only the traditional two-loop contributions with propagators and vertices carrying 
momenta p 1 ^ instead of p M , as depicted in the three graphs in fig 1. 

However, the linear terms in q in eqn 2.25 do contribute through a zero-momentum 
insertion, using the identity : 

- I X 

We have suppressed the C dependence in to simplify notation, and Q = Qo(0).This 
identity is of course nothing but the expansion of the 5 function in the defining equation 
2.4.* It gives us the zero-momentum insertions through the derivative. 

The identity 2.26 carries no volume factors, since the t has none right from the outset, 
see eqn 2.4 and above. Moreover, when it acts on a given term with a product of Q's, only 
terms linear in <5q(0) will survive the 5-function in 2.26. If there are n powers of q, then 
2.26 changes to one where n derivatives of Q act to the right. 

After these remarks we go through the usual procedure of expanding out of the 
exponential in 2.8 and doing the path integral, by contracting Q-fields. Apart from the 
terms, that give rise to the diagrams in fig 1, we have from the derivative in eqn 2.26 acting 

* The last step in this identity, the integration by parts, is strictly only permitted if the 
boundary conditions are such that the Q's are Gaussian variables. 
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on S( 3 ) and (S^) 2 : 



exp- 




s2 [^p^{^ )(c) ? Qo(p)Qo( - p) 



- 2t( 2 )(C)^Q (p)Qo(-p) 



2.27 



p 



V p 



v 





1 (^\c)) 2 J2Qo(p)Qo(-p) 



\\i^(C)\\ is the determinant of the matrix ^'(C) in the first term of eqn 2. 21. This deter- 
minant is the van der Monde determinant formed from the matrix C in eqn 2.16, as shown 
at the end of appendix D. 

Without knowledge of the explicit colour structure of eqn 2.27 the reader can easily 
check the following. All but one of the contractions in 2.27 are one loop contractions and 
hence of order 0(1). The only one that is 0(V) comes from the two-loop contraction we get 
from the second term in the straight brackets ; we have presented this term pictorially in 

(2) 

fig 2 and denote it by Up . The reader can easily verify this, using B.2 for the propagator. 

This contraction consists of two factors ; one is give by the renormalisation of the 
Polyakov loop : 




2.28(a) 



v 



The other is given by the zero-momentum insertion and equals : 



d_ 

dC 



U^\C) 



2.28(6) 



with U^\C) given in 2.23. Hence Uf> is the product of 2.28(a) and (b). 
Thus to two loop order U^ consists of two parts : 



u& = uf + 4 2) 
li 



2.29 



U^p is given by the graphs in fig 1 with the topology of the graphs contributing to the 
free energy but with the energies shifted through the background C, as described under- 
neath 2.21. U ( p ] is given by the 1-loop result 2.23 with the argument C shifted by the 
renormalisation of the loop 2.28(a) (see fig. 2). 

U^p is known to depend 6 ) on the gauge choice £. U p 2 ^ depends on the gauge choice 
only through 2.28(a), since 2.28(b) is independent of £ (see 2.23). Therefore there must 
be a relation between the gauge artifacts in both terms, in order to cancel and give a 
^-independent U^ 2 \ 

This relation is provided in the next subsection. 

2.e. BRST and gauge variation. 

In this subsection we will establish a relation between the BRST prediction for the 

(2) 

gauge - variance of Uj and the zero-momentum insertion 2.28 (fig 2(a)). In fact they are 
the same as a very simple argument will show. 

In fig 1 we have the three diagrams that contribute to U^ 2 \ In the £ gauge there 
are contributions linear, quadratic and cubic in (1 — £) from the propagators. The latter 
two are absent, as shown in appendix A. The terms linear in (1 — £) do not cancel. Rather, 
a simple observation shows they come in with a factor 3 from graph la, 2 from lb and 
1 from lc. Combining this with the combinatorial factors in front we see in fig 3(a) that 
U^p has a term linear in (1 — £) of the form : 

where II M „ is the one-loop self energy, and the momenta p lJ have colour shifted Matsubara 
frequencies defined underneath 2.21. 

The first sum is over the off-diagonal propagators, the second over the diagonal ones 
(in colour space). 

In the zero temperature case we obviously would obtain zero for the coefficient, because 
the BRST identities 10 ) tell us that the self energies are transverse. 

However, at non-zero temperature the BRST identities are not constrained by Lorentz 

(Euclidean) invariance. Although the BRST identities are local and therefore still valid at 

finite temperature their consequences are different. In our case, with a background field, 

they take the form 11 ) : 

n 6r 6r ' i i -i • 

= —— • — + ghost contribution + eventual matter contributions 

oQ oJ 

12 



r is the one-particle irreducible generating functional in an external field C ; at Q = it 
is related to Uf by V = f3VUf, since C is constant in space time. J is the external source 
coupling to the gauge variation of Q, the quantum field, through the term : 



J-D(Q + C)u 2.31 

The dot means summation over all degrees of freedom. 

To obtain information on the self energy IT in eqn 2.30 one derives the BRST identity 
with respect to the ghost-field uj and the quantum field Q ; setting all fields and sources 
to zero we get : 

n = . I . 2 32 

5Q5Q 5J5u 5Q 5J5u5Q 

This is the result to one-loop order. The first term in eqn 2.32 contains a double 
derivative with respect to the quantum fields. This double derivative equals the self en- 
ergy IT. All other terms are identically zero to this order. At zero temperature Lorentz 
invariance renders the last term in eqn 2.32 identically zero, and the one-loop self energy 
is transverse. However at finite temperature is not vanishing, and - by inspection- 

it equals the zero momentum insertion f3V 5 ^ } (see fig 3(c)) ; precisely through this zero- 
momentum insertion the renormalisation effect of the Polyakov loop did couple to the 
effective action, eqn 2.28. 

Note that only the temporal component survives in the second term of 2.32. 

So we have that both the gauge artifacts of the graphs in fig 1 and the full renor- 
malisation of the Polyakov loop couple to one and the same thing : the zero-momentum 
insertion of Qq into the 1-loop result. 

In order to make contact with the next section on the renormalisation of the Polyakov- 
loop we rewrite 2.32 in momentum space, with momenta and indices of the Cartan basis 
explicit : 

= ij#n& - ij^/^ 5 ^ l<i^j<N 2.33(a) 

and 

= ip^ d = 2,...,N 2.33(6) 

The colour-shifted momenta p lj are defined below 2.21. Qq(0) is the Euclidean space-time 
average of Qq(x) To obtain eqn 2.33(b) we used the fact that vanishes identically 

because of colour conservation. So only the off diagonal self-energies are non-transverse. 
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The terms linear in the gauge parameter (1 — £) in are proportional to the zero- 
momentum insertion in 2.33(a): 

3r =1>^ WM • 



When combining 2.33(a) and 2.34 we use the identity*) : 

N 

j?kl,lk,d ^ij,ji,d 



N 

£ f kl > lk > d f^> d = -(S ik + 8ji - 6u - S jk ). 



d=2 

(2) 

So finally substituting 2.33(a) into 2.30 we obtain for the coefficient of (1 — £) in Uj : 

Uf = Uf\t = 1) - 2(1 - ^Nj^B^C^B^Cij). 2.35 

Both in 2.34 and 2.35 we used eq. C.5 (with d=4 and k=l, 2 respectively) from 
Appendix C. 

The £ = 1 contribution in 2.35 is trivial to obtain and will be discussed in section 4. 

Let us note that the presence of fermions does not change the form of the BRST 
identity 2.32. So when fermions are present gauge artifacts can be treated the same way. 
The fermionic contribution shows up additively in the S3 factor in 2.35, as is clear from 
fig. 3 (c). 

2.f Comparison of gauge variation of and Up 

In this section we have found that both the renormalisation of the Polyakov loop, and 
the gauge artifacts in the free energy part of the effective action are coupled to the same 
object : the temporal and colour diagonal zero-momentum insertions into the one loop 
result 2.23 (see 2.28 (b) and the combination of 2.30 and 2.33(a)). The gauge variations 
of Up ' 1 and have the same structure. For the former we combine 2.28(a) and (b), for 
the latter 2.30(a) and 2.33(a). In section 3, eq. 3.18, we will see that evaluation of the 
renormalisation of the loop (2.28(a)) actually matches the gauge artifacts in 2.35, but with 
opposite sign. 

Therefore we will have complete cancellation of the gauge artifacts in the sum of the 
two, in two loop order. One expects this to continue for any order. 



*) As follows by inserting the definition D.2 for the structure constants and using the 
trace properties of the A lJ . 
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3. Renormalisation of the Polyakov-loop 



The renormalisation of the loop is not the renormalisation of the numerical value 
< ii >= y fdx < ti(x) >. 

Rather, what renormalizes is the relation between < t\ > and the particular choice 
of saddlepoint C that we made in eqn 2.16 in the phase of the loop. This renormalisation 
can - and will - contain a reflection of our gauge choice 2.12. 

As a general expectation one would say that the phase of the loop renormalizes with- 
out any extra ultraviolet infinities, beyond the usual ones that renormalise couplings and 
fermion masses. This should be so, because all these effects are thermal in nature. 

In fact it turns out, that the expectation values of the gauge invariant traces £i, fa, . . . tjv- 
indeed stay finite through one loop. From them one can reconstruct the finite renormaliza- 
tions of the eigenvalues of the loop easily. As a caveat let us look how proceeding through 
the renormalization of the gauge variant unitary matrix P(Aq) leads to a correction with 
unwanted properties.. 

So take the special unitary diagonal matrix P{Aq) = e l ^ c+9< ^° ,.It becomes after adding 
the one loop correction in fig 2. a a diagonal matrix which is not special, not unitary and 
not finite ! Nethertheless the expectation value of its trace t\ is finite as it should. The 
same holds true for the quantum average of the k th power of the Polyakov loop and taking 
its trace % (see appendix B). 

We have for the average of the loop matrix to order g 2 : 




dx < P(A ) > = e lCf3 - g 2 f dn [ 1 dr 2 e lCT2 

Jo Jo 



(Qo(T2)e^ ri - T2) Q (ri)e lC ^- ri )} + 0(/) 3.1 

where we used the definition 2.1 for P(Aq), 2.9 and 2.16. There is no x dependence after 
contraction of the Qo fields in 3.1, so we drop any reference to it for notational convenience. 

The order g 2 correction receives a contribution from all diagonal quantum fields Qq 
equal to : 

-9 2 ! dn r dr 2 e iC P (Q d (r2)Q d o(n)) 3.2 
Jo Jo 

because the diagonal Qq do commute with the diagonal e lCTi . Since 
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Aoo(T2-T 1 ) = (Qfe)g^(r 1 )> 

dp ( ! 



Po 



(2tt)3 V^+p 2 



(i-0 



(pl+P 2 ) 2 



3 «Po(t2-Ti) 



3.3 



we find with dimensional regularisation for the sum of all diagonal contributions to 3.2: 



2 1 N-1 (3 2 f d^p 
g 2 N 1 2 J (27T)-" 1 



4+2T V ^-A 00 (po,p) 



3.4 



The first term is zero ; the second as well, since Aoo is even in po. 

So we are left with the off-diagonal contributions Qq (t2)X ij ' (i ^ j). Let us denote the 
product of two A's by : 

X ij x ji = D u • _l ■ fixed 3 5 



and 



A"A" = />" % ^ j fixed 



3.6 



D" is the diagonal matrix, with all entries zero except the i th diagonal element, which 
equals 1/2. 

We use the identity 

e iCr X ij e -iCr = ^(Ci-C^ryj 3 _ 7 



and the propagator for the Qq excitation : 



,iP0(T2-Tl) 



3.8 



I 



Notice the occurrence of p$ = p + Ci — Cj. The time dependence of the propagator is of 
course periodic. 

With the help of 3.5 to 3.8 we can rewrite the 0(g 2 ) contribution in the form (see 
appendix B for details) : 



g 2 [{ (e-^-W - l] D u + (e« Ci - c 'W - l) D jj X A 



(2) 



i ' 



3.9 
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with 
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(27T)- 1 (p«)» 



The terms proportional D", D- 7 - 7 come from i < j , z > j respectively. 

The factors — in 3.10 stem from the integrations over the r- variables : they are 
there because of the non-locality of the loop. Observe in 3.9 that the term with r=l has a 
similar momentum structure as in 2.30, i.e. the loop on the left in fig 3.c proportional to 
Bi as in 2.35. We will see below, that the part with r=2 is projected out, after taking the 
trace of eq. (3.9). 

Obviously the renormalisation of the loop in 3.9 gives again a diagonal matrix and so 
eqn 3.1 can be rewritten as 

^ J dx< P(A )>=expi(C + g 2 5C)(3 + 0(g 4 ) 3.11 

with SC given by the matrix 

\ E f { (e- i{Ci - C ^ - l) D u + ( e *«*-W _ Aj 



l<i^j<N 



\2) 



+ (D u - />'-')A^ 



3.12 



We note that the trace of 5C is not zero ! It is : 



Tr5C = — ( cos (^ - C j)P ~ 1 ) A j J 2) 3 ' 13 

Note that is logarithmically divergent whereas is finite (see eqns C.5, C.15 and 
C.17). So our caveat at the beginning of this section has come true : as a matrix the 
expectation value of the Polyakov loop has become infinite, non-unitary and not special. 

One might argue, that the Polyakov loop as a matrix is not gauge invariant; in par- 
ticular under a periodic and constant (in x) gauge transformations 3.11 is similarity trans- 
formed, and so is the exponent C + g 2 5C. But the trace 3.13 stays invariant, so the infinity 
cannot be absorbed by the gauge transformation! 

On the contrary, when we first calculate the expectation value of the traces of the 
N — 1 powers of the loop, fi,?2 • • -fw-i ? there is no infinity, and the traces are real. To 
see this, let us start with SU(2) (N = 2). Then computing the trace of 3.9 gives us 

< t x >= cos/3(Ci - C 2 ) - g 2 \p^) sin(Ci - C 2 )(3 3.14(a) 
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with 

A;i) = ^(2 + (l-O)(^-0 3.14(6) 

The trace projects out the infinite terms ! The result 3.14 coincides with that of Belyaev 12 ) 
and that of ref 5. 

For SU(N) we find for the < ik > (see appendix B) : 

<i k >=^Tr((P(A )) k ) = ±Tr e* kC 

_ 9_pky T r f D ii_ D j3\ e <fcC A (v) fc = 1 . , iV - 1 3.15 

All infinite contributions do cancel in the trace. These N — 1 results in eqn 3.15 can be 
summarized by adding to C a diagonal traceless matrix 5C (not to be confused with the 
matrix 5C in 3.11 and 3.12), and to demand that the fe-th power of this matrix gives us 
the trace computed in 3.15 : 

< i k >= -^Trexpik(C + 5C)[3 + C(g 4 ) 3.16 

Identifying 3.16 with 3.15 gives for a diagonal element : 



3 

where B\{x) = x — \e{x). (See appendix C). 

From the anti-symmetry of in its argument CV, , it follows that the loop does not 
renormalize when C = 0, as expected. This is useful for knowledge of absolute minima 
(see section 6). It also shows that Yl!i=i^i = 0' i- e - the matrix SC is indeed traceless 
(and real and finite). 

(2) 

To find the full result for the Polyakov loop inserted Up we combine 2.28(a) and (b) 
with 3.17 to get: 

= 2(2+(l-£))g 2 Nj2Bi(C ij )B 3 (C ij ) 3.18 

i<j 

Eqns 3.17 and 3.18 are the main results of this section* . They generalize the results in 



* To obtain 3.18 more in line with the derivation of 2.35, take the insertion from 2.28 
as < ik >< a(f(i)fc^Qrf(o)) &int >• By inserting into 3.15 a complete set of diagonal 
generators one retrieves the matrix t^^ defined below 2.21. It factors out and therefore 
cancels the same matrix in the denominator of the average. Use of 2.34 and the subsequent 
identity for the f-symbols yields 3.18. 
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ref. 5 beyond the q- valley. They will be used when evaluating the full constrained effective 
action in the next two sections. 
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4. The gluonic constrained effective potential 



In this section we will give the one and two loop results. Everything is expressed in 
terms of the Bernouilli polynomials Bk(x), that we have compiled in appendix C. 
The result for an SU(N) gauge theory becomes to two loop order : 

U = U (1) + U {2) 4.1 



with 



and from 2.35 and 3.18: 



t/« = ^T 4 £s 4 (a,) 4.2 



i<j 



UW = uf(Z = l) + u¥\t = l), 4.3 

with the righthand side given in eq. (4.7) and (4.8). As anticipated all gauge artifacts 
have dropped out. The variable Cij equals 

C i:j = 4.4 
3 2tvT 

as in eqn 2.23, and B4, B2 are defined in eqns C.13 and C.14. 

A few remarks. LA 1 ) and are both of order A 2 for a generic value of Cij. However, 
for the specific choice in ref 5 : C\ = C2 = ■ ■ ■ = C^-i = 7^ and Cn = — Q we have 
only N — 1 terms : 

4tt 2 

Ud) = _T 4 (A - l)B 4 (q) 4.5 

U^ = ^.-5^L(N-l)B 4 ( q) 4.6 

So in this one parameter space the potentiel is O(N) smaller, and we will denote it 
by the " q- valley" . 

The reader can see this in the SU(3) example in fig 4. 

The q- valley formulae 4.5 and 4.6 were derived in ref 5 (eqn 5.14 in ref 5). Remarkable 
is the simplicity of 4.6: it just renormalizes multiplicatively the one loop result 4.5. This 
is lost outside of the q- valley. 

The actual derivation of in 4.3 is astonishingly simple in £ = 1 gauge. First we 
compute the contribution from the 3 graphs in fig 1. Using momentum conservation like 
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in A. 3, all integrands are of the form -Jp, where U , lj are two different momenta from 

i i 

the set of momenta Zi Z 2 and Z 3 flowing through the lines of fig la and b. The answer is : 

Uf = 9 2 Yl \f bca \ 2 {B 2 {C c )B 2 {Ca) - B 2 2 (0)) 4.7 



The colour indices 6, c, a run through the basis (ij) and d. C c = 2itT j if c = C c = 

if c = <i (<i labelling a diagonal generator). 

The contribution from the Polyakov loop inserion (£ = 1) is from 3.18 : 



U p 2) = Ag 2 N 



i<j 



4.8 



where the B^ are momentun integrals defined in appendix C and up to some multiplicative 
factors identical to the Bernoulli polynomials used before. Now the total Uf + Up reduces 
in the q-valley to 4.6, using the formulae in appendix C, and properties of the structure 
constants in appendix D. 

We would like to draw the readers attention to a remarkable, but unwanted property 
of Uf in 4.7. Let us write out its content in the q-valley: 



Uf = ^.^f^TWd - if - 29(1 - 4)] 4.9 



We see that the linear term causes an absolute minimum, where the value of U^ is 
negative. But remember the discussion in section 2.d: the free energy graphs alone follow 
from the 5 function constraint linear in the potential. As 4.7 is normalised by the the free 
energy graphs at C=0 we would expect the constrained path integral to give us a non- 
negative result! The fact, that it is not, just illustrates once more how gauge dependent 
effective potentials can give deceptive information. As another example 5 ) take this absolute 
minimum q m of U to two loop order. It will be of order 0(g 2 ). But by charge conjugation 
we get — q m and it is not hard to see, that there we have another absolute minimum: charge 
conjugation seems spontaneously broken! 

Both unwanted properties do disappear when we add 4.8, to get 4.6. More generally, 
also outside the q-valley our result 4.3 is non-negative and has only Z(N) minima. 
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5. Fermions 



In this section we add n / quarks (taken to be massless) in the fundamental represen- 
tation of SU(N). As fermions they are anti- periodic in time r. The action reads 

S = J dxdrq(0+4)q 5.1 

for every flavour. The Z(N) invariance is now broken through the boundary conditions: a 
gauge transformation which is periodic modulo a Z(N) phase will change the anti period- 
icity. 

We go through the same arguments as in section 2.b, so introduce background and 
fluctuation variables for the gauge fields. The next step is to compute the contribution 
Uq of the quarks to the constrained effective potential ; this has been done 1 ) and gives : 



5.2 



for rif quark flavours. The minus sign in front of 5.2 is due to the fermionic determinant, 
the Bernoulli function appears just as in the gluon case 4.2. Its argument is shifted over 
\ because of the anti-periodicity. The result is normalised to zero for Cj = 0. In what 
follows we will use the abbreviation : 

C{ = -^- + - 5.3 
2ixT 2 

To two loop order we will find-like we did for the pure gluon case in section 2.d - that 
the renormalisation of the Polyakov loop couples to the zero-momentum insertion of Uq . 
Likewise the arguments in section 2.e are valid for the linear gauge variation of the fermion 
case (see figsl, 2, and 3). 

The two loop result for n / quarks becomes : 



E { B 2 (C lJ )(B 2 (C()+B 2 (Cf)) - B 2 (C{)B 2 (Cf) 
+ ^ E( 2 ^(0)S 2 (C/ ) - Bl(C{)) 

i 

-|E s 3(^/)E^(^) + ( Ar2 - 1 ) 5 



144 



5.4 
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The terms in the braces come from the configuration in fig 1(d) where the gluon is 
off-diagonal. The next term stems from the diagonal gluons (see eqn D.7). The one but 
last term in 5.4 comes from the zero momentum insertion into 5.2 using eqn 3.17, and 
restores the minimum at Ci = 0. Thus, also in the fermion case charge conjugation and 
non-negativity are restored! But the Z(N) minima are no longer degenerate, because the 
symmetry is explicitely broken. 

In the g-valley both gluon contribution and fermion contribution are of order N, when 
q ~ N/2. This explains why the fermion contribution is comparable to the gluon contri- 
bution (fig 6). 

For 577(3) we plotted in fig 4 and 5 the full potential 4.2, 4.3 and 5.2, 5.4 with n/ = 0, 2 
respectively. The " q- valley" is seen on the border of the admitted values of the loop. It is 
plotted for n/ = 6 in fig 6. 
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6. Conclusion and outlook 



In this paper we established the perturbative expansion of the constrained effective 
action, to two-loop order in the large volume limit. Our result is consistent with the 
following prescription, to all orders : 

i) Compute the sum Uf(C) of all 1-PI free energy diagrams, with the energies shifted 
by C, and obtain 

U(C) = U f (C)-U f (0) 6.1 

ii) Compute the renormalisation of C to all orders to obtain 

C = C(C r ) 6.2 



iii) 

U(C(C r )) = U(C r ) 6.3 

is the final answer. 

Of course the result 6.3 is intuitively expected because of the relation 2.6 between U 
and the effective potential G 

As a function of C r , U will not depend on the gauge choice. U is therefore perfectly 
well defined in contrast to statements made by some authors. Of course, it may be that 
in three and higher loop order infrared divergencies will show up. They fall into two 
categories. The first one concerns divergencies due to the diagonal gluon propagators. 
The off diagonal propagators are protected by a value of CV, of order 1. The second 
category corresponds to small values of Cij of order g. The first category can be compared 
to the infrared divergencies in an Abelian L^l)^ -1 theory and is less severe. Some inroads 
into the latter have been made in ref 6). The second one is truly non Abelian and more 
severe. 

The other important issue was that the effective potential is non-negative and has 
Z(N) minima as absolute minima. A discussion of the latter can be found in ref. 13). Here 
we want to point out that the two are related. To this end, suppose 



U(C r ) > 
24 



6.4 



to all orders, because, at least formally, exp —(3VU(C r ) is a probability. We would like to 
know wether U(0) is an absolute minimum, i.e. wether U(0) = 0. According to 6.1, 6.2 
and 6.3: 



U(C r = 0) = 



6.5 



is equivalent to : 



F(C(C r = 0)) = F(0) 



6.6 



The simplest way to have 6.6 fulfilled is 



C(C r = 0) = 



6.7 



This is actually the case to two loop order (eqn 3.17). 

The authors of ref (13) have found a very ingenuous gauge choice, which they call 
Static Background Gauge (SBG) : 



In this gauge £ is fixed, £' — > 0. It looks very probable that in this gauge the Polyakov 
loop does not renormalize to any order of perturbation theory. If so, then indeed 6.5 is 
true to any order in perturbation theory.. 

With our gauge choice we found in section 3, that only traces of powers of loops are 
finite. The infinities in 3.12 are not ultra violet in nature ; they stem from the factors 
ij that originate in the r-integrations in 3.1, that is, in the non- locality of the loop. It 

(Po ) 

has consequences for the spontaneous breaking of charge conjugation as discussed in ref.5). 
Let us look first at the SU(2) case. Though the potential Aq becomes — A$ under charge 
conjugation (so in particular the colour diagonal matrix C will just flip sign), this will not 
affect the quantum average 3.14(a), since it is even in C. So the average of the trace of 
the loop is charge conjugation even, and no charge conjugation breaking effects can be 
measured with it. In general, for any N, charge conjugation leaves < P{Aq) > invariant 
as long as it is real. 

In between different vacua however the order parameter takes on complex values, 
(except for SU(2) of course), but this lies outside the scope of this paper. We will come 
back to the problem of computing the constrained effective action when the argument is a 
profile, as a function of one spatial variable 16 ). 




6.8 



25 



In fig 4 and 5 the reader will find the effective one loop potential for SU(S), with and 
without quarks. How the metastable state 15 ) in fig 6, near C = 4f , is sensitive to two loop 
contributions is shown in the same figure. 

The effects of lattice artifacts will be evaluated in ref 17). The Montecarlo evaluation 
could be done eventually with the multicanonical ensemble 18 ^ 19 ) ; this method avoids the 
difficulty of the peak between two vacua (see fig 6), which bars the MC-access to this 
domain of (/-values. Of course the massless fermions remain a practical problem. 

Some time ago 15 , it was pointed out that the total one loop effective action, i.e. the 
sum of the free energy at C=0, F(0), and U^\C) lead for high enough fermion number 
nf and appropriate values of C in the q- valley to absurd thermodynamical properties, in 
particular near the metastable minimum in fig 6. If we are only interested in local minima 
the presence of fermions is needed. 

But in a Montecarlo simulation one can apply an external field coupled to the Polyakov 
loop. 

In this latter context it is amusing to observe that the same happens for a system 
with gluons alone. For example, the maximum in fig 3 correponding to the center of the 
admitted values for SU(3), has negative entropy: explicitely, the free energy for iV 2 — 1 
gluons 1 ^ is — jgT7 2 T 4 (N 2 — 1). Adding to this the value of the one loop result for the effective 
potential 2.23 in the maximum gives for SU(2) f^T 4 and for SU(3) f^-2"" 4 , producing a 
negative entropy. This negative entropy may be a very serious default of the perturbative 
approach. Since at this order the free energy is just counting the degrees of freedom, 
the extremist would say the negative entropy indicates that we have overlooked hitherto 
unknown degrees of freedom in QCD. However, already the first derivative of the effective 
potential with respect to the order parameter shows peculiar behaviour: instead of the 
intuitively expected flat part- corresponding to the spontaneous breaking of the symmetry- 
it has smooth behaviour. This smooth behaviour is the opposite of what happens in Mean 
Field approximation, namely overcooling, and it fits with the fact that the second derivative 
of the effective potential is indeed convex (again, in contrast to Mean Field). Note that 
the effective potential is only convex in terms of the trace of the loop, not in terms of the 
background field C! Clearly some important physics remains to be understood. 
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Appendix A. Gauge terms proportional to (1 — £) 2 AND (1 — £) 3 



Gauge terms proportional to (1 — £) are treated in the main text, section 2.e. 

Terms proportional to (1 — £) 2 are involving contractions of the 3-vertices with two 
momenta. This contraction simplifies considerably the vertex. Consider diagram (a) in fig 
1. Its contribution to the (1 — £) 2 term equals : 



Y 2 I 2 {^pj Ci('a • hf - 2(/i • h){h ■ h)(h ■ h) + ll{h ■ hf) 

+ (1 -> 2, 2^3, 3 — > 1) + (1 — > 3, 2 -> 1, 3-^2)) A.l 



Because of the symmetry of the integration summation f ± 2 over li and I2 and colour, the 
two terms in A.l, not explicitely written, give the same result as the first one we will now 
further analyse. 

We use momentum conservation : 

h+h+h=0 A2 

to write 

2(k ■ lj) =l\- I 2 - I) A3 

for any triple (i,j,k) = (1, 2, 3). 

The first term / in A.l gives then : 

l= 9 — I (JL 1 A. — — - - — \ \ I 

48 J 1,2 V2^3 ^2 W'Z ^2 ^2. 



The last four terms add up to 



using 



9 —B 2 ■ B 2 A5 
48 

1 



/ 4 

1,2 1 



= A6 

2 



The first two terms give an equal contribution upon integration ; the second is seen to be, 
using A.2 : 

g 2 f (ll+2h-l 2 + l 2 ) _g* f (1 < 2h-h 1 , 



48 J lj2 ^1^2 48 7 1)2 Vt ^2 ^2 

The second term equals : 



g f 2/10/20 ^ g 



48 A, 2 ^2 
28 



because the correlation is odd in h (and Z 2 ). Using A. 4 to A. 8 we get : 

I=^(4B 1 -B 3 + 3B 2 -B 2 ) A.9 
The dots mean summation over colour degrees of freedom : 

Bi ' Bj = ^2 \fa,b,c\ 2 Bi(C a )Bj(C b ) 

a,b,c 

The third term in A.l is through exchange (1 — > 2) identical to the first one in A. 4 
we just computed. 

The middle term in A.l equals : 



/ (^3 + ^2 + ^1 ^1 ^2 h}\ ^3 
Jl,2 



t 

48 



"3 

The first term in A. 10 gives 



_ /4/2 _ /4,2 , 2/2/2/2N 1 A i n 



-rt 2 /" / 4 

^ ' /s All 



/ "I 

A,2 hi 



48 A,2 ^2 

whereas the other terms all cancel, using symmetry in the integration variables, and the 
result is for A.l : 4 

— ( 4Si • B 3 + 3B 2 • B 2 - - / A12 

The (1 — 2 contribution from diagram (6) in fig 1 can be worked out with the same tricks 
and gives the same result as in A.12, but with opposite sign. That is : the (1 — £) 2 term 
from fig 1 is zero. 

The (1 — £) 3 term is zero, because the 3- vertex vanishes when contracted with all three 
incoming momenta. 
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Appendix B. Quantum average of Polyakov loops 



In this appendix the trace of the quantum average of the n — th power of the Polyakov- 
loop is worked out. 

Let us introduce some notation. The matrix of order g 2 in eqn 3.1 is called L 2 : 

L 2 = -g 2 J P d Tl p dr 2 e* CT ^ J dxQ (r 2 x)e iC ^-^ 

Qo(Tix)e lC ^- T ^ B.l 

Remember C is a diagonal N x N matrix with i th entry Cj. 
Now we have for the propagator in a finite volume from 2.21 : 

(Qo{PoP)QoW)) = ^ P0 ,- P ' S P ,-py - (1 - ( ^f^ )4 ) B.2 



So : 



Po \ 



i / ^mnt)Qo(T 2 l)) -T^JJ^ - (1 - f)^, I e 

5.3 

in the infinite volume limit. 

In calculating Tr < L > we push e lCT2 through Q^ir^X 1 ^ to the right, which causes 
an extra phase e l ^ Ci ~ Cj ^ T2 to Qq (t 2 ). This phase adds to the Fourier coefficent e ipoT2 of 
Qofo) ; the same happens to the Fourier coefficient of Qq J (ti), it becomes e l ( p '°~( Ci ~ c ^ T2 . 

Having done this, integration over r 2 and t\ as in B.l, gives eqn 3.9 in the text. 

When computing the quantum average of TrP n (A), (n > 2) to order g 2 , we also need 
the first order matrix : 

L x = ig dTe tCT Q {T)e +lC ^- T) BA 
Jo 

To order g 2 we get (using P(A ) = e iCf3 + L x + L 2 ) : 

(TrP n (A)) =nTr (L 2 ) S n ~^ cp 

+ (n-l)Tr (LiLi) e l(n " 2)C/3 

+ (n-2)Tr(L 1 e* Cf3 L 1 )e^ n -V CfB 

+ ... + Tr(L 1 e i ^-^ c ' 3 L 1 ) + 0(g 4 ) B.5 
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We now push in each term the matrices e lCTi to the right, to have the two potentials next 
to one another. This leads to extra phases as below eqn B.l, and B.5 becomes, for a fixed 
pairiiig(gg'(75)Q5 < (7i)) = Ag ) : 



{TrP(Ar) =-g 



ip 13 (T2—Ti)iCij/3 inCP 



J dn dr 2 {nTrZT A* J e^ T2 " Tl V nC/3 

+ £ dr 2 dn { (n - l)Tr (d u A 1 ^ 

+ D 33 Aq J e ~ ipij ( T2 ~ Tl ^ iCi i P e inC P 
+ (n- 2)Tr tfi^i^-T^-nc^p jncp 

+D 33 A^ e~ ipij < - T2 ~ Tl )+ 2iC a $ e inC P j 



+ 



+D 33 AQ J e~ ipiJ ''(Tz-T^+iin-VdjPginCP^j 

Do the r-integrals and obtain, using the notation in 3.10 : 



B.6 



1 Tr(P n (A ))= <J 



N 



N 



(1- 
i 



T nTr (D u - D jj ) A\ 3 e inCp 

h 

+ n (TrD u ^ 2) e inCli (l - e~ lCi ^) 
+ TrD 33 Af 2) e* nCf3 (l - e* Ci ^)) 

+ Aj J 2) |l - e iCi ^\ 2 {(n - l)Tr (p^e™ e~ iC ^ + D jj e inC e iCi3 ') 
+ (n- 2)Tr (o"e inC e~ l2Clj P + D jj e inC e i2Ci3p ) 



+ 



B.7 



+Tr ^D ii e inC e~ i ^ n ~ 1 ^ Ci3 ' ^ + D 33 e inC ' e i ^ n ~ 1 ^ c%3 ^ 

The terms in the braces in B.7 do add up pairwise : the first and the last to n e t ^~ 1 ^ Cij/3 e' l ^^ Cx 
the second and the one but last to n e~*( %~ 1 ) Ci iP e l % {d+Cj)P ^ e t c ... Thus gathering all 
terms proportional to A( 2 )(u), we obtain as coefficient : 

n ^(d+c^p j e if _ e i(f-i)c i3 /3 + (2 - e iCi3 ^ - e ~ lCi313 ) 
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What is in the braces in B.8 adds up to zero. All what remains is the first term in B.7, 
which is nothing but eqn 3.15 in the main text. 
Amusingly, we can find another form of 3.17: 

This form is only true in special directions in the Lie algebra, given by the generalized 
hyper charges: 

Y k = —diag(k, k, ...k, k - N,k - N, - N). 

The entry k is integer, running from 1 to N-l. It is repeated N-k times, whereas k-N is 
repeated k times. In these directions the sign function has the wanted property 



J2«C il )-e(C jl )) = Ne(C ij ). 
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Appendix C Bernoulli polynomials 



The Bernoulli polynomials come about naturally in thermal field theory. We take 
from ref (14) 

^W-E^H^^cos^ C.l 

and 

2kB 2k _ 1 (x)=B' 2k (x). C.2 
Explicitely one has on the interval — 1 < x < 1 : 



B 4 (x)=x 2 (l-\x\) 2 -± 
3 1 

B 3 (x) = x 3 - -x 2 e(x) + -x C.3 
I 

B 2 (x) = x 2 — \x\ + — 



6 



Bx(x) = x- -e(x). 



e(x) is the sign function x/\x\. 

They are related to integrals of the type 



b -rVf d —l I r4 

d-2k- 2^ I ( 27r )d-i ((2nTn + 2nTx) 2 + P 2 ) k 



n 



no 

id-1 



R =T y / dd ~ lp (2nTn G + 2nTx) 

V d - 2k+1 -1 2^ I (27r)d _! {{2nTnQ + 2nTx) 2 + p2 )k 



B d = TJ2 [ 7^-^r(log((27rTn + 2nTx) 2 + P 2 ) - log((27rTn ) 2 + P 2 )) C.6 
no J ^ 



d is the number of dimensions of space time. For d = 4 we have 

2tt2 / _ ! 



B A {x) = ^-T 4 \B*(x) + ^ ) C.l 

B 3 (x) = —T 3 B 3 (x) C.8 

B 2 (x) = ^-B 2 (x) C.9 
33 



B x {x) = -^-B x {x) C.10 

4:7V 

Note that B4 and £?2 are even in x, B 3 and B\ odd. Hence £?2 and B\ are not analytic 
in x = 0. 

In deriving 4.3 from 4.7 and 4.8 we have used the identities : 

Si(x)S 3 (x) = B A {x) + ^B 2 (x) C.ll 



and 



(B 2 (x)) 2 = B 4 (x) + ^B 2 (x) + B 2 2 (0) C.12 



and 

B 2 (x) = B 2 (x) - 1 C.13 
Since it occurs frequently we define 

B 4 = B 4 + ^ C.U 

For the evaluation of A^ in eqn 3.10 it is useful to do the d-1 momentum integrations 
first, to find a ^-function of the type 14 ); 



00 



n=0 



(n + q) 2 



C.15 



with z = 4 — d. Using 14 ) 



SV ^ (z + l)(z + 2) z + 1 



valid for z non-negative, one finds 



T y f d P .. 1 — = J- Sl (c„) C.17 



n 



as in 3.14(b). 

The gauge term in A^ follows immediately from C.5 and C.10. 

ttm I I hm rr\ 4~ 1 \v> rv 1 

k (2) 



A/L will have terms like in C.15, but now with z=5-d. This gives us a pole 14 at d=4. 
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Appendix D Some group theory relations and two relations for the gluon poten- 
tial 



We introduce two unitarily related bases on the Lie algebra of SU(N). 

1) The Gell Mann basis A a , which is hermitean, traceless and normalised 

Tr\ a X b = -S ab D.l 

2) The Cartan basis X ij , X d (i ^ j), d = 2, 3, . . . , N, with : 

(X ij ) k i = -^SikSji , X d = -diag(l, 1, . . . , 1 - d, 0, 0, 0), r d = y/2d(d-l). 
V2 r d 

In both cases we define the structure constants / as : 

if abc = 2TrM a [M b ,M c ] D.2 

In basis 2) we have 

\fij,jk,ki\ = 2 

and 

N 

El/^| 2 = 1 DA 

d=2 

All other structure constants vanish in basis 2). 

It is useful to note that there are 2N(N — 1)(N — 2) /'s like in D.3, and that there 
are 3N(N — 1) sums like in D.4. 
So, in basis 1) we have 

J2\fa,b,c\ 2 = N(N 2 -l) D.5 

a,b,c 

and in basis 2) we find again : 

\f a ,b,c\ 2 = \ ■ 2N(N - 1)(N - 2) + 3N(N - 1) = N(N 2 - 1) D.6 



a,b,c 



as it should, because of the unitary relationship between the two sets. 

For the fermionic two-loop contribution 5.4 we need another identity for the coefficients 
r/s defined under D.l : 

*i \Ua 2 N 
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D.7 is true for all d = 2, . . . , N — 1. It is useful when computing the exchange of all diagonal 
gluons in fig 1(d). 

Let us finally relate the determinant H^^C)!! (appearing in the first term of eqn 
2.21) to the van der Monde determinant 

Jj (expiCi — expiCj) D.8 

l<i<j<N 

This equality is true up to a factor independent of C. This determinant was defined as 
the Jacobian between the variables ti(C),t2(C), ....tN-i(C), eqn 2.1(a), and the variables 
Cd = TrXdC, d=2,....N. From this it follows immediately, that a generic matrix element 
reads: 

t^(C) = ^L( eX pzA;Ci +expikC 2 + ... - (d - 1) exp ikC d ) D.9 

iV T d 

From this expression it is easy to see that the determinant of the (N-l)x(N-l) matrix 
tfc d (C) can be written as the determinant of the NxN van der Monde matrix, by adding 
and subtracting columns, and D.9 follows. 
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(a) (b) (c) (d) 

Fig 1: Graphs obtained from 2.12 with the topology of the free energy (a,b,c). 
The two loop fermion graph is shown in l.d. The symmetry factor is explicitly 
written. Continous (dashed) lines are gluons (ghosts). Constant background 
gauge is used, hence all energies are shifted through a constant amount. Their 

(2) 

contribution is given in 4.7 (for £ = 1) and 2.35 (if £ ^ 1) for the gluons, 
and in 5.4 for the fermions. 




Fig 2: Zero momentum insertion accompanied by renormalisation of the Polyakov- 

(2) 

loop, Up in eqn.2.28.The dot on the loops are the zero-momentum insertions 
2.28(b). 



(a) 



■ * o (o ■ o ■ < 

(b) (c) 

Fig 3: Relation between gauge variation of the free energy graphs (a), the gauge 
variation of the gluon self energy (b), and the zero-momentum insertion into 
the (cross-hatched) one loop free energy (c) through the BRST identity 2.32. 
The double bars represent the contraction of the momentum of the line into the 
vertex. The loop in the left part of fig. 3(c) equals the gauge variation of the 
renormalisation of the Polyakov loop in fig 2. 




0.8 




Fig 5: As in Fig 4, with rif = 2, using 5.2 and 5.4. 



1.0 2.0 3.0 

q 

Fig 6: The q- valley profile for the case N=3 and rif = 6. Solid curve is the 
one loop result 4.2 and 5.2., the dotted curve the one and two loop result at 
a s =0.1. The dashed curve is the pure gluon result. 



